Lubrication effects on the flow of wet granular materials 
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We investigate the dynamics of a partially saturated grain-liquid mixture with a rotating drum 
apparatus. The drum is partially filled with the mixture and then rotated about its horizontal 
axis. We focus on the continous avalanching regime and measure the impact of volume fraction 
and viscosity of the liquid on the dynamic surface angle. The inclination angle of the surface is 
observed to increase sharply to a peak and then decrease as a function of liquid volume fraction. 
The height of the peak is observed to increase with rotation rate. For higher liquid volume fractions, 
the inclination angle of the surface can decrease with viscosity before increasing. The viscosity where 
the minima occurs decreases with the rotation rate of the drum. Limited measurements of the flow 
depth were made, and these were observed to show only fractional changes with volume fraction 
and rotation speeds. We show that the qualitative features of our observations can be understood 
by analyzing the effect of lubrication forces on the timescale over which particles come in contact. 

PACS numbers: 45.70.Ht, 45.70.Mg 



I. INTRODUCTION 



Granular flows have been studied over the last few 
decades in order to confront the challenges posed by sev- 
eral unresolved questions ^, 'S, . Most of these studies 
have focused on the relatively simple dry granular system 
with air as the interstitial fluid. However, lately, there 
has been a keen interest in understanding the dynamics 
of granular flows wetted by a small amount of liquid con- 
sidering their prevalence in industry and nature 0, Q. 
The addition of a small amount of liquid to a dry granu- 
lar system results in the formation of liquid bridges which 
binds the particles. The static behavior of these liquid 
bridges has been studied for quite some time [1, 0| on 
the grain scale and the cohesive force generated is shown 
to be dependent on the surface tension of the liquid and 
grain size. On the macroscopic scale, there are a few 
quantitative descriptions of the impact of the liquid on 
the overall behavior. A wet pile of granular material is 
more stable than a dry pile [8, 9] and the failure in a wet 
pile occurs in the bulk rather than at the surface as for 
the dry case [l^, [HI- Ratholing and jamming can oc- 
cur during hopper drainage under conditions where dry 
grains would normally flow. Furthermore, particle segre- 
gation can be su ppr essed or enhanced by making them 
cohesive [H, [Tsl Il4 [T5| . Bond number (Bo) which is 
given by the ratio of the strength of the capillary bridge 
to the weight of a grain has been introduced as a useful 
measure of the cohesivity of the material [l^ . 

An important system to study granular matter is a 
partially filled horizontally rotated drum. This system 
has been extensively used to study stability, avalanching, 
mi xing and segregation properties of dry materials [ITI . 
m, [S US mi, m m. The flow is mostly confined to 
the surface and is only a few layers deep. A transition 
from stick-slip to continuous avalanching is observed as 
the rotation rate is increased. More recently, the rotating 
drum apparatus is being used to investigate wet granular 
materials because it offers many important features [ill . 



lid I24 l25l . |26| . First, because the system is periodic, 
it can be investigated over long times in a steady state. 
Second, by simply varying the rotation rate, various time- 
scales can be accessed. An attractive practical feature is 
that the system can be sealed very easily, which prevents 
contamination and evaporation of the liquids. 

Recently, using this system, some of the static prop- 
erties of a wet pile has been well characterized. Mea- 
surements of the angle of maximum stability as a func- 
tion of the volume fraction of the liquid, the grain size, 
and the dimensions of the wet pile inside the drum have 
been made [26[. Two different models have been devel- 
oped to describe the observed dependence. The first phc- 
nomenological approach considers the internal friction of 
the grains and models the increased stability due to in- 
crease in cohesive force [13, HI]. The second approach 
considers the geometric stability of the grains and the 
increases stability due to the orientation of the capillary 
bridges [2^. Both approaches yield similar scalings for 
the maximum stability in terms of the liquid and grain 
properties which are observed to be in good agreement 
with experimental observations [26j . 

Tegzes et al [2^ l25j have performed detailed measure- 
ments of wet granular avalanches in the stick slip regime. 
They found that the nature of the avalanches depended 
on the volume fraction of the added liquid. The three 
flow regimes: granular, correlated, and viscoplastic were 
characterized by the nature of the stick-slip avalanches. 
Recently, simulations have considered the effect of adding 
simple cohesive forces on granular flow [2^ [s^] ■ 

In this paper, we discuss our investigation of the effect 
of liquids in the continous avalanche regime. We measure 
the angle of inclination as a function of rotation rate of 
the drum, and the volume fraction and viscosity of the 
added liquid. A wide drum is used to minimize the effect 
of the side walls. We define a lubrication number which is 
given by the ratio of the time needed for the two particles 
in neighbouring streamlines to move towards each other 
to the time needed to pass each other without contact. 



2 




FIG. 1: An image of the experimental apparatus. The drum 
is partially filled with a mixture of glass beads and a small 
amount of silicone oil and rotated about its horizontal axis. 
The dotted box shows the region captured in the image to 
measure the surface angles. 



This lubrication number is found to be useful in describ- 
ing the observed variations of the inclination angles. 



II. EXPERIMENTAL APPARATUS 

The experimental apparatus consists of a clear plexi- 
glass drum of diameter 28.5 cm and width W = 14.5 cm 
rotated about its horizontal axis as shown in Fig. [TJ The 
drum is rotated using a computer controlled stepper mo- 
tor allowing for a rotation rate (w) variation from 0.00456 
rpm to 13.69 rpm. The drum is 40% filled with spher- 
ical glass beads of diameter d = 1 mm, premixed with 
small amount of viscous liquid. The amount of liquid 
added is reported in terms of volume fraction Vf , defined 
as the ratio of the volume of liquid to the volume of dry 
glass beads. The volume of dry beads is calculated as the 
weight of the beads divided by the bead density {pb = 2.5 
g cm~'^). The liquid used in the experiments is Silicone 
oil with four different viscosities (v — 0.05,0.2,1.0 and 
10.0 cm^s-i). 

Images are taken from one end of the drum and in 
the central region as indicated by a dotted box in Fig. [1] 
(The image resolution of 1024 x 768 pixels corresponds to 
a region of 18 cm x 13.5 cm). Images are captured using 
back lighting to identify the flowing surface. We have 
found that such a lighting scheme is less sensitive to a 
thin layer of grains which stick to side walls of the drum 
and therefore capture the surface of the granular mixture 
away from the side walls. The region with grains appear 
black in the image and unfilled regions of the drum ap- 
pear white due to the back light. The free surface pro- 
file is extracted using an edge detection algorithm from 
which the angle of the surface is obtained. Measurements 
are made for the steady fiow achieved after initial few rev- 
olutions of the drum. The reported values for the surface 



angle are obtained by averaging over 500 images (several 
revolutions of the drum). During the measurements, the 
drum is rotating continuously and particles are always 
in motion so that the particle-liquid mixture remains ho- 
mogeneous throughout. The surface fiow of the parti- 
cles exhibits a stick-slip motion at lower rotational rates 
which changes to continuous motion at higher rotational 
rates. In the remainder of the paper, we discuss the flow 
behavior in the continuous avalanching regime. 

III. RESULTS AND DISCUSSION 

We first discuss the variation in the shape of the flow- 
ing surface with changing rotational speeds and the de- 
pendence of surface angles properties on the choice of 
location along the flowing surface where the measure- 
ments are made. The effect of interstitial liquid volume 
fraction and viscosity on the surface inclination is quali- 
tatively discussed next. We then examine the relative ef- 
fect of the liquid induced forces for various experimental 
parameters in order to explain the observed fiow behav- 
ior. 



A. Inclination angle of the flowing surface 

The surface of the pile is almost linear at low to and be- 
comes progressively non-linear for higher w. Figure [21^ a) 
shows a typical surface profile at low and high cu for wet 
systems. The non-linearity of the surface at higher uj re- 
sults in different slopes along the surface. Figure [2jb) 
shows the slopes of the surface in the uppermost and 
lowermost portion of the extracted surface profile. The 
slopes (surface angles) are obtained by averaging over a 
region in which the profile is approximately linear. The 
surface angle of the upper section increases with uj while 
that in the lower section decreases slightly with uj. This 
behavior is observed for dry case as well and the wet case 
for all the viscosities and volume fractions of the liquids 
investigated. Within the range of experimental param- 
eters studied, the higher surface angle corresponding to 
the upper section appears to be the most sensitive to the 
variation in system parameters and in further discussion 
we focus on this angle to examine the effect of volume 
fraction and viscosity of the liquid added to particles. 

B. Effect of liquid volume fraction and viscosity 

Figure [3] shows the surface angle (6) plotted against 
the rotational speed for different volume fractions (Vf) 
of the added liquid of viscosity ly — 1.0 cm^s~^. Also 
shown is the data for the dry granular system (filled sym- 
bols) for the entire range of rotational rates studied. For 
all the cases studied, the surface angle increases with 
UJ, but the rate of increase progressively reduces with in- 
creasing amounts of liquid added. Notice particularly the 
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FIG. 3: (a) Variation of surface angle (6) with the rotation 
rate for different volume fractions (V/) of the added liquid 
of viscosity u = 1.0 cm^s~^. Filled circles denote the corre- 
sponding data for dry system. 



FIG. 2: (a) Free surface profile of the flowing particles as ob- 
tained from the images for two different rotational rates, (b) 
The surface angles calculated from upper and lower sections of 
the the free surface as a function of rotation rate. The data is 
obtained for the wet system with the added liquid of viscosity 
{u = 1.0 cm^s"^) and volume fraction {Vf = 0.48 x 10"^). 

nearly flattened profile at the highest Vf . This rate of in- 
crease with rotational speed for all the volume fractions 
is slower than the quadratic dependence observed for the 
dry systems ■ The added liquid tends to bind the par- 
ticles strongly resulting in lower surface angles and slower 
rate of increase with rotational speeds. A small amount 
of liquid added {Vf = 0.08 x 10~^) shows a marginal, 
though a curious change. The angles for the dry sys- 
tem are higher than the wet system at lower rotational 
rate with the trend reversing at higher rotational rates. 
We attribute this behavior to the slight fluctuations in 
the angles for dry system due to humidity effects which 
were not controlled in the experiments. Further increase 
in the volume fractions results into higher surface an- 
gles than those for the dry case, suggesting the influence 
of the acting capillary attractive forces between parti- 
cles thereby increasing the resistance to flow. However, 
adding more liquid also increases the lubrication between 
particles causing them to flow much easily as is seen for 
the highest volume fraction and higher rotational rates, 
where the angles observed are actually lower than those 
for the corresponding dry case. Similar behavior is ob- 
served for liquids with different viscosities. 

To gain a further understanding about the effect of 
added liquid, we replot the data as shown in Figure 3] 
for four different viscosities studied. Consider flrst the 
case for ly = 1.0 cm^s~^. For lower rotational speeds, 9 
shows a rapid increase with volume fraction followed by 
a steady increase over the remaining range. The increase 
in the angle with volume fraction is a direct consequence 



of the cohesive capillary force due to the liquid bridges 
formed between particles as mentioned above. An impor- 
tant feature which is revealed is the formation of a peak 
in the surface angle proflles at higher rotational speeds. 
The angle rapidly increases to a maximum and then de- 
creases continuously to a much lower value at higher vol- 
ume fractions which is due to the flow being lubricated 
by the interstitial liquid. Similar trends are seen for dif- 
ferent viscosities as well. The peak, however, is broad for 
lower viscosities and narrower for higher viscosities. A 
minimum observed at intermediate volume fractions and 
the highest rotational rate for i/ = 10.0 cm^s^^ is not 
quite clear and suggests a complex interplay of the forces 
due to the liquid. 

To further elucidate the effect of viscosity, we plot the 
surface angles versus the viscosity for different rotational 
speeds and for the highest volume fraction studied as 
shown in Figure [5] For lower rotational speeds the an- 
gles increase monotonically with viscosity which is ex- 
pected considering that the resistance to flow would in- 
crease with viscosity. However, the angles show a reverse 
behavior at the highest rotational rates, wherein the liq- 
uid actually assists the flow and reduces the angles. A 
combination of these effects is seen at intermediate rota- 
tional rates wherein the angles decrease to a minumum 
and then increase for higher viscosities. 

The discussion so far suggests that the effect of the 
added liquid is sensitive to the rotational rates employed. 
The primary effect of adding the liquid is to generate 
two newer forces between the particles viz. the attrac- 
tive capillary forces and the repulsive viscous forces in 
addition to the gravity and fricitional forces which exist 
in a dry system. In the next section, we examine these 
forces and their relative dominance with respect to each 
other for different rotational rates which would help in 
understanding the observed behavior of the surface an- 
gles. 
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FIG. 4: Variation of the surface angle with volume fraction for 
different rotation rates, (a) v — 0.05 cm^s~^, (b) (a) v — 0.2 
cm^s"\ (c) V = 1.0 cm^s"^ and (d) u = 10.0 cm^s"^ In all 
cases, the rotational rates range from 1.37 rpm to 13.69 rpm 
as specified in (d). 



45 



40 



35 



30 



J (B = 2.28 rpm 
^ CO = 4.56 rpm 
(B = 6.85 rpm 
(B = 9.13 rpm 
11.41 rpm 
13.69 rpm 




0.1 1.0 

V (cm7s) 



10.0 



FIG. 5: Variation of the surface angle with viscosity for dif- 
ferent rotational rates and fixed volume fraction (V/ = 6.0 
xlO"^). 



C. Effect of liquid induced forces 

When a small amount of liquid is mixed with dry 
beads, the liquid wets the bead system and forms liquid 
bridges between the particles bonding them together by 




FIG. 6: Schematic of a liquid bridge formation between two 
spherical particles. 9 is the liquid/solid contact angle, is 
the half-filling angle, R =(d/2) is the diameter of the parti- 
cle and h is the distance between two grains, ri and r2 are 
respectively the radius of the neck and the meridian profile. 



a weak attractive capillary force. A schematic of such a 
liquid bridge between two particles is shown in Figure [S] 
In the limit of small amount of liquid, the gravity effects 
can be neglected and the total attractive force can be ex- 
pressed as a sum of the component of surface tension (cr) 
force along the length of the bridge and the force aris- 
ing due to the pressure difference between air and liquid. 
Further, assuming d/2 ^ r2 ^ ri and h <C 2ri cos 6* and 
that the liquid bridge is approximately cylindrical (flat 
profile), a simplified expression can be obtained relating 
the total force to the bridge volume {V) and separation 
distance {h) which is given as 31 1 



Fcap ~ Trdcr cos 9 



4V 



(1) 



where V = nd[H^ib) - h'^]/A, H{r) = h + 2r^/d, r is 
the radial distance from the center of the bridge and b is 
the radius of the wetted area. The capillary force varies 
directly as the volume of the liquid bridge which depends 
on the volume fraction of added liquid. 

The pressure change in the liquid bridge arises due to 
the relative displacements of two particles with respect 
to each other while in flow and the two can be related by 
the Reynolds equation which can be simplified to deter- 
mine the expression for the viscous force acting on the 
particles. For finite volumes of liquid (cylindrical bridge) 
the viscous force is given as [3l| 



X 1,. 



-Tvpivd 



1 



1 dh 

h~di' 



(2) 



where pi is the density of the liquid, dh/dt is the rate of 
approach of particles towards each other. H{h) is calcu- 
lated from the expression given above. The viscous force 
varies directly as the rate of displacement of particles, is 
negligible for large separation distances and diverges as 
particles come very close to each other. 

Direct measurments or estimates of these liquid in- 
duced forces in addition to the gravity and frictional 
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FIG. 7: An image taken with 1/4 second exposure shows the 
depth of the surface flow. Volume fraction Vf = 0.08 x 10~^, 
rotation rate w is 2.28 rpm. The arrows indicate the general 
direction of the flowing particles. 



forces for various experimental parameters would eluci- 
date the observed behavior of surface angles. However, 
our experimental techniqe prevents the direct estimates 
of these forces between the particles. So instead, we ex- 
amine different scenarios determining the relative signif- 
icance of different forces with repect to each other. This 
can be done by calculating the appropriate timescales 
corresponding to the motion of particles relative to each 
other which can be obtained through velocity profiles in 
the flowing layer of particles. The spreading of the inter- 
stitial liquid on the end plates combined with the sticking 
of the particles to the walls prevents accurate measure- 
ments of the velocities in the flowing layer. To get around 
this, we instead measure the flowing layer thickness (6) 
using streakline photography similar to that used ear- 
lier [2l| for the dry systems. This involves taking time ex- 
posed photographs to obtain particle streaklines as shown 
in Fig. [T] Because we are seeking estimates, we focus on 
the central flow region and calculate the maximum layer 
thickness over there. We note here that the problems 
faced while measuring accurate particle velocities do not 
hinder the measurements of the flowing layer thickness. 
The typical flow direction of particles and the location 
for layer thickness (S) measurement is shown in Fig. [T] 

The variation of the flowing layer thicknes with rota- 
tional rate (w) for different volume fractions and for one 
particular viscosity of the liquid is shown in Fig. [SI The 
layer thickness changes very slowly with rotational rate, 
which is similar to the results for the dry systems with 
the measurements using NMR [H, and particles vi- 
sualization ^2l\ techniques. The rate of change decreases 
with increasing rotational rates, particularly for higher 
volume fractions with the layer thickness approaching a 
constant value at high enough rotational rates, similar to 
the trends observed for dry systems. Further, the layer 
thickness increases much slowly with rotational rates for 
higher volume fractions and interestingly the layer thick- 
ness at the highest volume fraction becomes almost con- 
stant and independent of the rotatational speed. Similar 
behavior is also observed for liquids with different viscosi- 
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FIG. 8: Variation of the normalized flow depth (S/d) with 
rotation rate for different volume fractions of the added liquid 
{u = 0.2 cm^/s), where d is the diameter of the particle. Filled 
circles denote the corresponding data for dry system. 



ties. For the steady flow of particles in the continuous 
flow regime, the flux Q of the particles (at the central 
location where S is measured) can be written from mass 
conservation as 



(3) 



where pt is the density of glass beads (assumed to be ap- 
proximately equal in the flowing layer and rotating bed) , 
Rq is the radius of the drum, z is the length of the cylin- 
der and Vx is the average velocity in the flowing layer. 
The average velocity in the flowing layer can be then 
written as 



LuRi 



26 



(4) 



Using Eq. 21 we can thus get the estimates of the repre- 
sentative velocities in the flowing layer for various rota- 
tion rates, volume fractions and viscosities of the inter- 
stitial liquid. For the particles flowing in the x-direction 
(see Figure [7]), we denote t.^ — d/{Vx) as the time re- 
quired for the two particles in the neighbouring stream- 
lines to travel past each other, where d is the particle 
diamete r. Substituting for Vx from Eq.|4] gives 



2dd 



(5) 



The timescale tx primarily varies inversely as the rota- 
tional rate (cj) since 6 is only weakly dependent on uj and 
volume fractions as seen from Figure [S] 

We now define ty as the timescale for the same two 
particles to approach each other while squeezing out the 
liquid from between. Balancing the viscous force between 
the two particles given by Eq. [5]to the weight of the par- 
ticle {ird'^ p(,g/6) and then integrating the resulting equa- 
tion gives 



_ 9 vpi 

^11 — 



4 dgpb 



2b^d{h, - ho) 



{2b^ + h,d){2b^ + hod) 



6 



In 



1 + 2b'^/hsd 
l + 2b^/hod 



(6) 



where hg is the surface roughness of the beads, and for 
the glass beads is approximately 10~'^d. b is the radius 
of the wetted area of the particle. The above equation 
gives an estimate of the time required by the two particles 
separated by a distance hg to approach each other to a 
distance within surface roughness (hg). 

The ratio Ra = tx/ty, after substituting for the exper- 
imental parameters, can be then written as: 



Ra 



0.13 S/{luv) 



(2b'^ + hsd)(2b'^+bd) 



In 



l+2b'2/h,d 
l+2b/d 



(7) 



where, we have assumed that /iq ^ i-e. liquid bridges 
just before or after snapoff are semi-spherical droplets on 
the surface of the beads. We can now investigate various 
scenarios under which different forces can be dominant. 
For Ra ^ 1, the two particles pass each other without 
coming in contact (thus minimizing the frictional forces) 
while for Ra ^ 1 the viscous forces between the particles 
dominate the overall flow behavior. In the following we 
now estimate these time scales for different rotational 
rates, volume fractions and viscosities of the liquids used 
in our experiments. 

For small Vf, as when the beads first get uniformly 
coated with the liquid, and liquid bridges are formed be- 
tween the grains, 6 ~ /is [231 ■ In that limit, both the 
terms in the denominator of Eq. ([7]) tend to and there- 
fore Ra 1 and particles touch as they flow past each 
other and kinetic friction/ viscous forces become impor- 
tant resulting in the higher surface angles. This flow 
regime corresponds to the surface angles increase with 
volume fraction (see Fig. [4]) leading to a maximum. The 
peak in the profile would in that case correspond to the 
maximum value of Ra. 

On the other hand, for large enough volume fractions, 
b ~ d/2, ho d/2 and using hg and S ~ IM 

(from Fig. [5]), we get Ra 10^'^ / (tov) . For highest 
rotational rate (w = 13.69 rpm) and lowest viscosity 
{f = 0.05 cm2 s-1), Ra O{10-^) < 1. We thus ex- 
pect the surface angles to be lower and is indeed seen 
from the decrease in surface angles with volume frac- 
tion (see Fig. [D^a)). This effect is more pronounced at 
the highest viscosity used {ly = 10.0 cm^ s~^) for which 
Ra — > 0(10^'^) <^ 1. Note the faster decay in the surface 
angles with increasing volume fractions in Fig. [4Kc,d). 
The rotation rate and viscosity, thus, act together to give 
maximal lubrication. This is reflected in the non-intuitive 
decrease in the surface angles with increasing viscosity at 
highest rotational rate (see the topmost curve in Fig. [5]). 

The simple analysis presented here based on the time 
scales for the particles to come in contact serves to ex- 
plain most of the effects due to varying volume fractions 
of the liquid for different rotational rates and viscosities. 
The primary dependence on the viscosity of the liquid or 
the rotational rates for a given volume fraction, however. 



cannot be explained using this analysis. For example, the 
angle increases with viscosity at lower rotational rates, 
the reason being the increase in the drag forces with vis- 
cosity leading to additional resistance to flow and hence 
higher angles (see the bottom curve in Fig. [5]). Using 
Eq. ([7]) alone in this instance would cause the correspond- 
ing Ra values to change from 0(1) to O(10~^) suggest- 
ing a decrease in the angle contrary to that observed. A 
similar argument also holds true for the case where the 
angles increase with rotational rates (the corresponding 
Ra actually decreases in those instances). Apart from 
these behaviors, the analysis cannot as well account for 
the more complex behaviors like the sharp minimum in 
the surface angles observed at higher rotational rates and 
highest viscosity (see Fig. HJi) or the decrease in the angle 
followed by an increase with increasing viscosity observed 
for intermediate rotational rates (see Fig. [5]). Under- 
standing these complicated effects would require a more 
rigorous study involving estimation of the different com- 
peting forces which is beyond the scope of this work. 



IV. SUMMARY 



In summary, we have investigated the effect of volume 
fraction, viscosity of the liquid, and rotation rate of the 
drum on the inclination angle of a glass bead pile in a 
rotating drum. The surface angle show a non-linear in- 
crease with the rotational rate for any particular liquid 
volume fraction and viscosity. For any rotation rate and 
liquid viscosity, adding small amounts of liquid increases 
the surface angle than the corresponding values for the 
dry system which is a direct consequence of the liquid 
bridge formation between particles thus increasing the 
resistance to flow. However, continuing the addition of 
liquid eventually leads to a constant slow increase in the 
surface angles at low rotational rates whereas the angles 
show a surprising decrease at high enough rotation rates 
even though the strength of liquid bridge does not de- 
crease. The effect is seen to be more dramatic for higher 
viscosities. By doing an analysis of the time scales related 
to the approach of two particles towards each other dur- 
ing flow using we derived a criterion for determining the 
relative effects of friction and lubrication between flowing 
particles. The analysis is valid primarily for accounting 
the effects due to varying liquid volume fractions and 
captures most of the observed variation in the angle of 
inclination. However, we are unable to make simple com- 
parisions of the scaling because of different competing 
forces in the experiments which are not measurable using 
the techniquies employed here. Additionally, it would be 
helpful to have detailed velocity profiles and flow depths 
away from the sidewalls as a function of liquid content, 
and will be the focus of future work. 
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